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The calculation of electronic conductance of nano-scale junctions from first principles is a long 
standing problem in molecular electronics. Here we demonstrate excellent agreement with exper- 
iments for the transport properties of the gold/alkanediamine benchmark system when electron- 
electron interactions are described using the many-body GW approximation. The main difference 
from standard density functional theory (DFT) calculations is a significant reduction of the contact 
conductance, Gc, due an improved alignment of the molecular energy levels with the metal Fermi 
energy. The molecular orbitals involved in the tunneling process comprise states delocalized over 
the carbon backbone and states localized on the amine end groups. We find that dynamical screen- 
ing effects renormalize the two types of states in qualitatively different ways when the molecule 
is inserted in the junction. Consequently, the GW transport results cannot be mimicked by DFT 
calculations employing a simple scissors operator. 



I. INTRODUCTION 

The conductance of a molecule sandwiched between 
metallic electrodes is sensitive to the chemical and elec- 
tronic structure of the molecule as well as the detailed 
atomic structure of the metal- molecule contact. Varia- 
tions in the contact geometry beyond experimental con- 
trol lead to an undesired spread in measured conduc- 
tance properties. For the most commonly used thiol an- 
choring group, these effects are rather pronounced due to 
the many possible contact geometries resulting from the 
large strength of the Au-S interaction. Amine groups 
have been shown to produce more well defined trans- 
port properties^ which can be understood from the rel- 
atively weak AU-NH2 bond leading to larger structural 
selectivitjJ2]. 

Even for a given junction geometry, a quantitatively 
accurate description of electron transport from first prin- 
ciples remains a formidable task. Numerous studies 
based on density functional theory (DFT) have shown 
a significant overestimation of conductance relative to 
experimental value^^'^ (an exception to this trend oc- 
curs for small molecules, like H2^^ and CO^^ which 
chemisorp strongly to the electrodes resulting in resonant 
transport through broad, partially filled resonances). 
The inability of DFT to describe off-resonant tunnel- 
ing in the simplest molecular junctions limits the pre- 
dictive power of the DFT-based approach to qualitative 
trends. It is now broadly accepted that the failure of 
DFT is mainly due to its wrong description of the molec- 
ular energy levels. Indeed, physically motivated correc- 
tion schemes have shown that much improved agreement 
with experiments can b e obt ained after shifting the DFT 
molecular energy level^^'^. Such corrections are sup- 
posed to remove the self-interaction errors inherent in 
standard DFT exchange-correlation (xc) functions and 
account for image charge effects induced by the metal 
contacts. The drawback of the approach is that is as- 
sumes a weak coupling between molecular orbitals and 



metal states and treats the image plane position as a free 
parameter. 

The (self-consistent) GW approximatiorfi^l^ which is 
rooted in many-body perturbation theory, was re- 
cently found to yield a considerable improvement 
over DFT for the conductance of gold/benzenediamine 
junction^SSl Physically, the GW approximation corre- 
sponds to Hartree-Fock theory with the bare Coulomb 
interaction v = l/|r — r'| replaced by a dynamically 
screened Coulomb interaction ly(aj) = e^^(w)u. In con- 
trast to standard DFT, the GW approximation is self- 
interaction free and includes image charge effects due to 
the metal contacts through the correlation part of the 
self-energy^i2al^ As a consequence, it provides quantita- 
tively accurate predictions of energy gaps in systems with 
highly di verse screening properties ranging from isolated 
molecule^2iIIIl over semiconductors^^, to metals^^. This 
property becomes particularly important for a metal- 
molecule interface where the electronic structure changes 
from insulating to metallic over a few Angstrom. 

In this work we use the GW approximation to study 
the role of exchange-correlation effects for the energy 
level alignment and electron transport in short alkane 
chains coupled to gold electrodes via amine linker groups. 
The gold/alkane junction is a benchmark system for 
molecular charge transport and h ave been exhaustively 
investigated experimentalljP^^'^lHlIl_ focus here 

on the amine-linked alkanes to avoid the uncertain- 
ties related to the gold-thiol contact geometry which is 
presently under debat^^^HiII_ We note that very recently 
it was shown that alkanes can be bound directly to gold 
electrodes without the use of anchoring groupdrSSl. 

The transport mechanism in (short) saturated molec- 
ular wires is coherent tunneling via molecular orbitals 
with energy far from the Fermi energy. The conductance 
versus chain length {N) thus follows an exponential law 
of the form 

G = G,exp(-/3JV) (1) 
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Recent experimentally reported values for the decay 
constant /3 of alkanediami ne /gol d junctions are in the 
range 0.9-1.0 per C atorrpJi^lMJ^ but earlier measure- 
ments also showed values around 0.^^. Although pre- 
vious studies based on DFT have yielded /3 values 
within the experimental range, the contact conductance, 
Gc, is typically overestimated by around an order of 
magnitude^ "^'^^ By comparing DFT and GW calcula- 
tions for 7V-alkanediamine molecules with iV = 2, 4, 6 we 
show that the wrong Gc values are a result of incorrect 
level alignment in the DFT calculations. Indeed, GW 
yields a Gc in close agreement with experimental values. 
We find a pronounced orbital and length dependence of 
the quasiparticle (QP) corrections to the DFT energies 
resulting from the different shape and localization of the 
molecular orbitals. The QP corrections range from —0.5 
to —2.5 eV and can be qualitatively explained from a 
classical image charge model. 

II. METHOD 

The junction geometries were optimized using the 
real space projector augmented wave method GPAWSma 
with a grid spacing of 0.2 A and the PBE functional for 
exchange and correlation (xc)^^. The molecules were at- 
tached to Au(lll) surfaces, modeled by an eight layer 
thick 4x4 slab via small four atom tips as shown in 
Figure [l] The surface Brillouin zone was sampled on a 
4x4 Monkhorst pack fc-point grid, and the structures 
including molecule, Au tips, and outermost Au surface 
layers were relaxed until the residual force was below 
0.03 eV/A. We have considered iV-alkanediamine junc- 
tions with = 2, 4 and 6. Key structural parameters can 
be found her^SSI_ Por calculations of the molecules in the 
gas-phase we include 16 A of vacuum between molecules 
in repeated super cells. 




FIG. 1: Supercell used to model the gold/alkanediamine junc- 
tions. Similar supercells have been used for N = 2 and 
4 (not shown) and key structural parameters can be found 
hereP. The black box indicates the extended molecule region 
within which the GW self-energy is evaluated. The inclusion 
of the Au tips in the extended molecule is sufficient to include 
screening effects from the electrodes on the molecular energy 
levels. 

All transport calculations where performed using the 
method described in detail in Ref.f^. In brief we employ 
a basis set of numerical atomic orbitals corresponding to 
double-zeta plus polarization (DZP) for the Au atoms 



and double-zeta (DZ) for the atoms of the molecules. 
We use rather diffuse basis functions with a confinement 
energy shift of 0.01 eV. This ensures that the calculated 
work function of Au(lll) and the Kohn-Sham energy lev- 
els of the molecular junction are within 0.1 eV of those 
obtained from accurate grid calculation^2o!^ xhe trans- 
mission function is obta ined from the Meir-Wingreen 
transmission formulgP^^ 

T(e) = Tr[G^{e)TL{e)G''ie)Tn{e)]. (2) 

The retarded Green function of the extended molecule is 
calculated from 

G'^iE) = [{E + ty)S -Ho + V^c - AVii[G] 

~j:L{E)-^RiE)-i:,c[G]{E)]-' (3) 

Here 5, Hq, and Vxc are the overlap matrix, Kohn-Sham 
Hamiltonian and the local xc-potential in the atomic or- 
bital basis, respectively. 77 is a positive infinitesimal. 

The lead self-energy, Si/i? , incorporates the coupling 
to the left and right electrodes and is obtained using stan- 
dard techniques--. The term AVh gives the change in 
the Hartree potential relative to the DFT Hartree poten- 
tial already contained in Hq. Finally, the last term is 
the many-body xc self-energy which in this work is ei- 
ther the bare exchange potential, Vx, corresponding to 
Hartree-Fock or the GW self-energy. As indicated both 
the Hartree potential and the xc self-energy depend on 
the Green function. The latter is evaluated fully self- 
consistently using a simple linear mixing of the Green 
functions. The GW self-energy is evaluated for the ex- 
tended molecule (indicated by the box in Figure [ij. How- 
ever, only the part corresponding to the molecule is used 
while the remaining part is replaced by the DFT xc- 
potential. This is done to include non-local correlation 
(image charge) effects from the electrodes in the GW 
self-energy of the molecule while preserving a consistent 
description of all metal atoms at the DFT level. We have 
verified that the calculations are converged with respect 
to the size of the extended molecules, see Ref.l^for more 
details. We represent all energy dependent quantities in 
Equation ([3]) on a large energy grid ranging from —200 eV 
to 200 eV with an energy grid spacing of 0.01 eV. 



III. RESULTS AND DISCUSSION 

A. Energy level alignment 

The alignment of the molecular energy levels relative 
to the electrode Fermi level is of great importance for the 
transport properties of molecular junctions and seems to 
be the dominating effect at low bias voltage. At higher 
bias voltages many body calculations on small model sys- 
tems suggest that electron correlations induce additional 
shifting and broadening of the molecular levels which can 
also affect the transport properties^-. Here we focus on 
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TABLE I: Calculated HOMO and HOMO-2 energies aligned 
to the vacuum level and in units of eV. 



Method 


Orbital 


i V — z 


AT = 4 iV = 6 


DFT-PBE HOMO 


-4.9 


-5.1 


-5.1 




HOMO-2 


-8.5 


-8.2 


-8.0 


HF 


HOMO 


-10.2 


-10.5 


-10.5 




HOMO-2 


-13.3 


-12.9 


-12.8 


GW 


HOMO 


-8.5 


-8.6 


-8.6 




HOMO-2 


-12.2 


-11.8 


-11.6 



TABLE II: Calculated HOMO and HOMO-2 energies in the 
junction relative to the electrode Fermi level. 

Method Orbital N = 2 N = 4 N ^6 



DFT-PBE HOMO 


-4.3 


-4.2 


-4.4 




HOMO-2 


-6.5 


-5.9 


-5.7 


HF 


HOMO 


-8.1 


-8.1 


-8.3 




HOMO-2 


-10.4 


-9.9 


-9.7 


GW 


HOMO 


-4.8 


-4.9 


-5.2 




HOMO-2 


-8.4 


-8.2 


-8.2 



the low bias regime and postpone the finite bias effects 
to a later study. 

The molecular orbitals (MOs) of the alkanediamine 
chains comprise states which are delocalized over the car- 
bon backbone and states which are localized on the NH2 
end group. We shall consider the highest occupied molec- 
ular orbital (HOMO) and HOMO-2 as representatives 
for the two classes of states, see Figure [2]j a). We note 
that the HOMO-1 is similar to the HOMO with slightly 
lower energy given by the coupling of the two end groups 
across the wire. In Table |T] and |lT] we list the energy 
of the HOMO and HOMO-2 calculated with DFT-PBE, 
Hartree-Fock (HF) , and GW for the molecules in the gas- 
phase and junction, respectively. In the gas-phase, all 
three methods predict the HOMO energy to be almost 
independent of molecular length. This is clearly due to 
its end group localized character. In contrast, the en- 
ergy of the HOMO-2 level shifts upward in energy as 
the molecular length increases. This reflects its extended 
nature and can be interpreted as a band discretization 
effect. To the best of our knowledge no experimental 
results exists for the ionization potential of alkanedi- 
amine molecules. However, for the closely related butane 
molecule {N = 2 alkane with CH3 end groups) we ob- 
tain a GW calculated HOMO energy of —11.4 eV in very 
good agreement with the experimental ionization poten- 
tial of 11.2 eVfSD. In comparison, the DFT-PBE HOMO 
energy is severely underestimated at —7.9 eV. This find- 
ing agrees well with pr evious studies on a larger range of 
small moleculed20|24i25j 

In the junction, the molecular orbitals, \ipn), have been 
obtained by diagonalizing the DFT Hamiltonian corre- 
sponding to the molecule. The projected density of states 



(PDOS) of such a state is then given by the spectral func- 
tion, — l/7rIm(?/'„|G"'(i?)|V'n), where G is the appropri- 
ate Green function (calculated with DFT, HF, or GW). 
The level position is defined as the first moment of the 
PDOS. Figure [2](b) shows the PDOS for the HOMO and 
HOMO-2 for the 6- alkanediamine junction as calculated 
with DFT-PBE (upper panel) and GW (middle panel). 
The lower panel shows the PDOS obtained from a DFT 
calculation where the molecular levels have been shifted 
to match the GW levels, i.e. after adding to the Kohn- 
Sham Hamiltonian a generalized scissors operator of the 
form, SgfT ^ ^j^QP _ ,DFT)|^^^^^^|, Here e^^ de- 

note the QP energy obtained from the GW calculation. 
We see that the main features of the GW spectral func- 
tion can be well reproduced by the shifted DFT Hamil- 
tonian although small differences remain. A similar con- 
clusion was reached in ReLl^for a gold/benzenediamine 
junction. The molecular orbital energies from a GW cal- 
culation include the dynamical response of the electron 
system to an added or removed electron through the cor- 
relation part of the self-energy. In general correlations 
tend to shift the filled levels up and the empty levels 
down relative to the bare Hartree-Fock energies. This 
is because the inclusion of screening reduces the energy 
cost of removing /adding electrons to the molecule. When 
a molecule is brought into contact with a metallic junc- 
tion its environment changes from insulating to metal- 
lic. This implies extra screening of an added or removed 
electron which will cause the filled levels to shift upwards 
and the empty levels to shift downwards even more than 
for the isolated molecules, i.e the gap will shrink rela- 
tive to its gas-phase value. It has been shown previously 
that DFT in (semi)local approximations and Hartree- 
Fock completely misses this important effeclP3t2ll, 

In Figure ^c) and (d) we show the QP corrections to 
the HF and DFT Kohn-Sham energy levels as function 
of molecular length. The results for the HOMO and the 
HOMO-2 are denoted by circles and squares. We first 
notice that the QP corrections are very significant with 
absolute values reaching almost 4 eV with a pronounced 
orbital and length dependence. The Hartree-Fock QP 
corrections are all positive showing that HF places the 
occupied levels lower than predicted by GW. This is in 
contrast to the corrections to the DFT levels which are 
all negative in agreement with the well known underes- 
timation of ionization potentials by DFT. In contrast to 
Hartree-Fock the Kohn-Sham QP corrections are smaller 
for molecules in the junction rather than in the gas-phase. 
In fact the DFT HOMO level position is relatively close 
the GW level position and only lie 0.5-0.8 eV higher. The 
fact that the DFT-PBE description of molecular energy 
levels is much better in the juncti on than in the gas- 
phase agrees with previous finding^^SHHlll^ a,nd can be 
explained from the origin of the PBE functional in the 
homogenous electron gas^-. 

It is instructive to consider the shift in the molecular 
energy levels due to correlation effects coming from the 
metal electrodes. In simple terms this corresponds to the 
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FIG. 2: (a) Isosurfaces for the HOMO and HOMO-2 orbitals 
of the 6-alkanediamine molecule, (b) HOMO and HOMO-2 
PDOS in the junction obtained from DFT-PBE (top), GW 
(middle) and DFT with a generalized scissors operator (bot- 
tom), (c) Quasiparticle corrections to the Hartree-Fock levels 
in the gas-phase (open symbols) and in the junction (filled 
symbols) as a function of molecular length A''. The HOMO 
and HOMO-2 are denoted by circles and squares, respectively, 
(d) Same as (c) but for DFT rather than Hartree-Fock. 



shift induced by image charge effects. In order to isolate 
the part of the correlation energy originating from the 
metallic electrodes we define the quantity, 



^junction 



(^n )gas-phasc (4) 



-HF 



which is shown in Figure [3| The result can be under- 
stood qualitatively by considering a classical model to 
account for the screening effect of the electrodes. Clas- 
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FIG. 3; Change in the correlation energy of the HOMO and 
HOMO-2 energy levels when the molecules are taken from the 
gas-phase into the junction. This value represents the shift 
of molecular levels due to the enhanced screening provided 
by the metallic electrodes (an image charge effect in simple 
terms). 



sically a charge distribution between two perfect met- 
als will experience an image potential. For example, the 
image potential for a point charge halfway between two 
metal surfaces separated by a distance L is ~ 10.0/L (eV 
A)^^. This predicts a oc 1/L dependence for the image 
charge effect. The Au tip atoms in our simulations are 
about 8 A apart for the 2-alkadiamine junction giving 
a rough estimate of 1.3 eV for the image charge effect 
in qualitative agreement with the GW calculations. The 
HOMO experiences a larger image charge effect than the 
HOMO-2 which can be understood from the fact that its 
charge density is located closer to the metallic surfaces. 
In the limit of an infinitely long wire the HOMO-2 will be 
spread out over the entire molecule and the image charge 
effect should vanish. On the other hand in this limit the 
HOMO will stay localized near the surface and therefore 
approach a non-zero constant image charge potential. If 
we model the HOMO charge density as a point charge of 
half an electron on each of the amine groups we can esti- 
mate this limiting constant to be 3.6/(2c?) (eV A), where 
d is the distance to the nearest metal surface. Taking 
d to be around the Au-N bond length (2.34 A) gives a 
limiting value estimate of 0.8 eV. Again, this seems to be 
in qualitative agreement with our GW findings. 

Finally, we discuss the coverage dependence of the en- 
ergy level position for alkanediamine-Au junctions. It 
was shown in Ref.^^ that the DFT level position for amine 
linked molecules is strongly dependent on coverage. In 
contrast to the screening (image charge) effects discussed 
above, which appear in the correlation part of the self- 
energy, this is a purely electrostatic effect resulting from 
the localized surface dipoles formed at the AU-NH2 bond. 
To investigate the dependence of the energy levels on cov- 
erage for our junctions we have performed DFT calcula- 
tions for a range of transverse supercell dimensions for 
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FIG. 4: The molecular valence band edge (or HOMO) as 
a function of coverage for an A'^ = 6 alkane. The numbers 
indicate the number of surface atoms. 



the geometry shown in Figure [T] In Figure |4] we show 
the PDOS of a methylene unit in the central part of the 
molecule for transverse super cells with 2 x 2 to 8 x 8 sur- 
face atoms. The PDOS peaks and band edge is shifting 
up in energy as the transverse supercell size is increased 
in agreement with the results of Ref.^®. The inset shows 
the energy shift obtained from reading off the shift in the 
PDOS as a function of coverage defined as 77 = 1/A^surfacc, 
where iVgurfacc is the number of surface atoms. For a su- 
percell dimension of size 3x3 and larger, the shift is seen 
to be directly proportional with the coverage as expected 
for a two dimensional array of dipoleJ^. This allow us 
to to extract the electrostatic shift corresponding to the 
single molecule limit. We find that the electrostatic en- 
ergy shift when going from a 4 x 4 supercell to the single 
molecule limit is indeed significant with a value around 
1 eV. 



B. Transport calculations 

The transmission function of the 2, 4 -and 6- 
alkanediamine junction geometries were calculated using 
the GW and the PBE xc potential as approximations for 
Sxc in Eq. ([3]). To include the coverage dependence we 
have simulated the low coverage limit 77 = 0"*" by perform- 
ing calculations using the 4x4 junction (corresponding 
to ?7 = 1/16) with all molecular levels shifted up by 1 eV 
using a simple scissors operator self-energy. 

The transmission function calculated using GW is 
shown in Figure [5] on a logarithmic scale. The transmis- 
sion functions for different molecular lengths have very 
similar shapes in the important region near the Fermi 
level Ep^ however, the magnitude is increasingly sup- 
pressed as function of molecular length. The similarity 
of the transmission functions may at first seem surprising 
since we have shown that the position of the molecular 
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FIG. 5: The transmission function calculated using GW for a 
molecular length of N=2, 4 and 6. 
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FIG. 6: Calculated conductance plotted as a function of the 
molecular length for coverages corresponding to 4 x 4 Au 
atoms per molecule (77 = 1/16) and extrapolated to the single- 
molecule limit (77 = 0). 



energy levels show some length dependence. In partic- 
ular the HOMO level was found to decrease in energy 
by 0.5 eV when N increases from 2 to 6 (see Table [TT| . 
This shift is indeed visible in transmission function in the 
range -4.0 to -6.0 eV where the HOMO is located. On 
the other hand the feautures in the transmission func- 
tion around the Fermi level are determined by the local 
electronic structure of the Au tips . 

The zero bias conductance is obtained from the trans- 
mission function at the Fermi level, G — GoT{Ep) where 
Go = 2e^/h is the unit of quantum conductance. The 
zero bias conductance is plotted in Figure [6] as function 
of molecular length. We have also included the DFT re- 
sults for comparison. The dashed lines show the best fits 
to the exponential form Gat = Gc exp(— /3iV). The values 
for f3 and Gc corresponding to the single-molecule limit 
(77 = 0) and 4x4 Au atoms per molecule (7/ = 1/16) 
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DFT 




GW 


Exp. 


V 


1/16 





1/16 





Gc 


0.11 


0.10 


0.02 


0.04 0.02°, 0.035'= 


P 


0.99 


0.92 


1.01 


1.02 0.97", 0.93', 0.91" 



"Ref.EH 
=Ref.El 

TABLE III: Calculated contact conductance Gc in units of 
Go and exponential decay constant P per carbon atom. 77 
denotes the coverage. 



are reported in Table III together with experimental val- 
ues. The rather weak effect of coverage on the conduc- 
tion properties is in agreement with the findings in Ref.''^ 
where a 4-alkanediamine Au junction with 3x3 and 4x4 
surfaces was considered. Our DFT results are in rea- 
sonable agreement with previous DFT studies showing 
decay factors in the range 0.83 — 1.01 and conductance 
resistances in the range (0.09 - Q.2S)G^^. While the 
/3 values obtained with GW are rather close to the DFT 
calculated ones, the contact conductance is reduced by a 
factor of 3 — 5 depending on the coverage. This is a di- 
rect result of the molecular levels lying further away from 
Ep (by 0.5-2.5 eV, see Figure [2jd)) in GW compared to 
DFT. 



correlation effects for the energy level alignment and low- 
bias conductance of gold/alkanediamine molecular junc- 
tions. Based on many-body GW calculations we found 
that the origin of the overestimation of the contact con- 
ductance, Gc, by standard DFT is due to wrong energy 
level alignment in the junction. The absence of self- 
interaction and the inclusion of image charge screening ef- 
fect through the GW self-energy improves the desription 
of the energy levels and yield values for Gc and the de- 
cay constant j3 in good agreement with experiments. The 
quasiparticle corrections to the DFT energy levels showed 
a significant orbital dependence ranging from —0.5 eV 
to —2.5 eV due to the different shape and localization 
of the molecular orbitals. Our results demonstrate that 
quantitatively accurate calculations of conductance from 
first-principles are feasible, although computationally de- 
manding. 
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